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We investigate in detail the qualitative behaviour of the class of Bianchi type B spatially ho- 
mogeneous cosmological models in which the matter content is composed of two non-interacting 

Q\ . components; the first component is described by a barotropic fluid having a gamma-law equation of 

^\ ' state, whilst the second is a non-interacting scalar field <j> with an exponential potential V(<j>) = Ae fe ^ . 

ON , In particular, we study the asymptotic properties of the models both at early and late times, paying 

particular attention on whether the models isotropize (and inflate) to the future, and we discuss the 

i— ( ' genericity of the cosmological scaling solutions. 
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,—l ' I. INTRODUCTION 

CO ' 

iy~j | Scalar field cosmology is of importance in the study of the early Universe and particularly in the investigation 
of inflation (during which the universe undergoes a period of accelerated expansion Jl|,|§l). One particular class of 
inflationary cosmological models are those with a scalar field and an exponential potential of the form V(4>) = Ae k ^, 
where A and k are non-negative constants. Models with an exponential scalar field potential arise naturally in 
alternative theories of gravity, such as, for example, scalar-tensor theories. 

Scalar-tensor theories of gravitation, in which gravity is mediated by a long-range scalar field in addition to the 
usual tensor fields present in Einstein's theory, are perhaps the most natural alternatives to general relativity (GR). 
In the simplest Brans-Dicke theory of gravity (BDT; ||), a scalar field, <f>, with a constant coupling parameter u)q, 
5h i acts as the source for the gravitational coupling. More general scalar-tensor theories have a non-constant parameter, 
^P" w(<p), and a non-zero self-interaction scalar potential, V{4>). Observational limits on the present value of ujq need not 
^ , constrain the value of u> at early times in more general scalar-tensor theories (than BDT). Hence, more recently there 
has been greater focus on the early Universe predictions of scalar-tensor theories ofgravity, with particular emphasis 
on cosmological models in which the scalar field acts as a source for inflation [§Jj]. BDT ( and other theories of 
gravity, such as, for example, more general scalar-tensor theories and quadratic Lagrangian theories and also theories 
undergoing dimensional reduction to an effective four-dimensional theory [5]), are known to be conformally equivalent 
to general relativity plus a scalar field having exponential-like potentials 5 M ■ 

Scalar-tensor theory gravity is currently of particular interest since such theories occur as the low-energy limit 
in supergravity theories from string theory and other higher-dimensional gravity theories j7j. Lacking a full non- 
perturbative formulation which allows a description of the early Universe close to the Planck time, it is necessary 
to study classical cosmology prior to the GUT epoch by utilizing the low-energy effective action induced by string 
theory. To lowest order in the inverse string tension the tree-level effective action in four-dimensions for the massless 
fields includes the non-minimally coupled graviton, the scalar dilaton and an antisymmetric rank-two tensor, hence 
generalizing GR (which is presumably a valid description at late, post-GUT, epochs) by including other massless fields; 
hence the massless bosonic sector of (heterotic) string theory reduces generically to a four-dimensional scalar-tensor 
theory of gravity. As a result, BDT includes the dilaton-graviton sector of the string effective action as a special case 
(u> = — 1) 0. String cosmology has recently been investigated by various authors ||, and, in particular, Q presented 
a qualitative analysis for spatially flat, isotropic and homogeneous cosmologies derived from the string effective action 
when a cosmological constant term is included. A discussion of how exponential potentials arise in effective four- 
dimensional theories (in the so-called conformal Einstein frame) after dimensional reduction from higher-dimensional 
theories such as string theory and M-theory is given in JTo| | . 

A number of authors have studied scalar field cosmological models with an exponential potential within GR. 
Homogeneous and isotropic Friedmann-Robertson- Walker (FRW) models were studied by Halliwell || using phase- 
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plane methods (see also || ). Homogeneous but anisotropic models of Bianchi types I and III (and Kantowski-Sachs 
models) have been studied by Burd and Barrow ||llf in which they found exact solutions and discussed their stability. 
Lidsey [jiJJ and Aguirregabiria et al. found exact solutions for Bianchi type I models, and in the latter paper a 
qualitative analysis of these models was also presented. Bianchi models of types III and VI were studied by Feinstein 
and Ibanez [ fL4| , in which exact solutions were found. A qualitative analysis of Bianchi models with k 2 < 2, including 
standard matter satisfying various energy conditions, was completed by Kitada and Maeda |fl5| . They found that 
the well-known power-law inflationary solution is an attractor for all initially expanding Bianchi models (except a 
subclass of the Bianchi type IX models which will recollapse). 

The governing differential equations in spatially homogeneous Bianchi cosmologies containing a scalar field with an 
exponential potential exhibit a symmetry |16| , and when appropriate expansion- normalized variables are defined, the 
governing equations reduce to a dynamical system, which was studied qualitatively in detail in jl7| |. In particular, the 
question of whether the spatially homogeneous models inflate and/or isotropize, thereby determining the applicability 
of the so-called cosmic no-hair conjecture in homogeneous scalar field cosmologies with an exponential potential, was 
addressed. The relevance of the exact solutions (of Bianchi types III and VI) found by Feinstein and Ibanez p4| , 
which neither inflate nor isotropize, was also considered. In a follow up paper 18 the isotropization of the Bianchi 
VII/! cosmological models possessing a scalar field with an exponential potential was further investigated; in the case 
k 2 > 2, it was shown that there is an open set of initial conditions in the set of anisotropic Bianchi VII^ initial data 
such that the corresponding cosmological models isotropize asymptotically. Hence, scalar field spatially homogeneous 
cosmological models having an exponential potential with k 2 > 2 can isotropize to the future. However, in the case 
of the Bianchi type IX models having an exponential potential with k 2 > 2 the result is different in that typically 
expanding Bianchi type IX models do not isotropize to the future; the analysis of |ll| indicates that if k 2 > 2, then 
the model recollapses. 

Recently cosmological models which contain both a perfect fluid description of matter and a scalar field with an 
exponential potential have come under heavy analysis. One of the exact solutions found for these models has the 
property that the energy density due to the scalar field is proportional to the energy density of the perfect fluid, 
hence these models have been labelled scaling cosmologies With the discovery of these scaling solutions, it has 

become imperative to study spatially homogeneous Bianchi cosmologies containing a scalar field with an exponential 
potential and an additional matter field consisting of a barotropic perfect fluid. The scaling solutions studied in 
|^,|lj, which are spatially flat isotropic models in which the scalar field energy density tracks that of the perfect 
fluid, are of particular physical interest. For example, in these models a significant fraction of the current energy 
density of the Universe may be contained in the scalar field whose dynamical effects mimic cold dark matter. 

In 1 22 1 the stability of these cosmological scaling solutions within the class of spatially homogeneous cosmological 
models with a perfect fluid subject to the equation of state p 7 = (7— l)p 7 (where 7 is a constant satisfying < 7 < 2) 
and a scalar field with an exponential potential was studied. It is known that the scaling solutions are late-time 
attractors (i.e., stable) in the subclass of flat isotropic models [ p0|j2l] . In |^2| it was found that that the scaling 
solutions are stable (to shear and curvature perturbations) in generic anisotropic Bianchi models when 7 < 2/3. 
However, when 7 > 2/3, and particularly for realistic matter with 7 > 1, the scaling solutions are unstable; essentially 
they are unstable to curvature perturbations, although they are stable to shear perturbations. Although these solutions 
are unstable, since they correspond to equilibrium points of the governing dynamical system, the Universe model can 
spend an arbitrarily long time near these scaling solutions, and hence they may still be of physical importance. 

In addition to the scaling solutions described above, curvature scaling solutions and anisotropic scaling solutions 
are also possible. In |^3| homogeneous and isotropic spacetimes with non-zero spatial curvature were studied in detail 
and three possible asymptotic future attractors in an ever-expanding universe were found. In addition to the zero- 
curvature power-law inflationary solution and the zero-curvature scaling solution alluded to above, there is a solution 
with negative spatial curvature where the scalar field energy density remains proportional to the curvature, which 
also acts as a possible future asymptotic attractor. In [^3J spatially homogeneous models with a perfect fluid and a 
scalar field with an exponential potential were also studied and the existence of anisotropic scaling solutions was also 
discovered; the stability of these anisotropic scaling solutions within a particular class of Bianchi type models was 
discussed. 

The purpose of this paper is to comprehensively study the qualitative properties of spatially homogeneous models 
with a barotropic fluid and a non-interacting scalar field with an exponential potential in the class of Bianchi type 
B models (except for the exceptional case Bianchi Vl^/g), using the Hewitt and Wainwright formalism p^ , Pq| . In 
particular, we shall study the generality of the scaling solutions. The paper is organized as follows. In section II we 
define the governing equations, which are modified from those developed in |2E) |, and discuss the invariant sets and 
the existence of monotonic functions. In section III, we classify and list all of the equilibrium points, and their local 
stability is discussed in section IV. We give a detailed analysis, including hctcroclinic orbits, for a subset of Bianchi 
type Vlft, models in section V. We leave conclusions and discussion for section VI. 
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II. THE EQUATIONS 



We shall assume that the matter content is composed of two non-interacting components. The first component is 
a separately conserved barotropic fluid with a gamma-law equation of state, i.e., p = (7 — where 7 is a constant 
with < 7 < 2, while the second is a noninteracting scalar field 4> with an exponential potential V(<p) = Ae fe ^, 
where A and k are positive constants (we use units in which 8ttG = c = 1). By non-interacting we mean that the 
energy-momentum of the two matter components will be separately conserved. 

The state of any Bianchi type B model with the above matter content can be described by the evolution of the 
variables 



H, <t+,(t, 6, a, n + ,4>, cpj E R , (2.1) 

where the evolution of the state variables are given as equations (5.8) and (7.8) in Wainwright and Ellis |^5[ with the 
addition of the Klein-Gordon equation for the scalar field, 

ip + iH^j) + kV{(t)) = 0. (2.2) 

By introducing dimcnsionless variables, the evolution equation for H decouples and the resulting reduced system 
has one less dimension pj|. Defining p5|,fLT| 



£ -f± Z--L A-— A- — 



n + = ^±, v = ^-, $ = ^M, n = -| 5l (2.3) 

+ H' VGH' V3H 3ff 2 ' V ' 
the differential equations for the quantities 

X = (£+,£, A, i,JV+,tf,$) e M 7 (2.4) 

are as follows: 

S' + = ( <Z -2)S+-2iV, (2.5) 

£' = 2(g-2)£-4AiV + -4£ +J 4, (2.6) 

A' = 2(g + £+ - 1) A + 2(£-iV)iV + , (2.7) 

A' = 2(g + 2£ + )i, (2.8) 

= (g + 2£+) iV+ + 6 A, (2.9) 

*' = (g-2)#- i\/6fc$ 2 , (2.10) 

$' = (g + l + -V / 6fc*)$ ) (2.11) 

where a prime denotes differentiation with respect to the time r, where dt/dr = H. The deceleration parameter q is 
defined by q = — (1 + H' /H), and both iV (a curvature term) and (a matter term) are obtained from first integrals: 

(7 = 2S^ + 2i] + i(37-2)f7 + 2* 2 -$ 2 , (2.12) 

N=\n%-\\A\ (2.13) 

Q = 1 - * 2 - $ 2 - £^ - ±-N- A. (2.14) 
The evolution of Q is given by the auxiliary equation 

n' = 51(2(7-37 + 2). (2.15) 

The parameter I = 1/h where h is the group parameter is equivalent to Wainwright 's h in p5| ]. If I < and A > 
then the model is of Bianchi type VT/j. If / > and A > and iV + 7^ then the model is of Bianchi type VII/j. If 
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I = then the model is either Bianchi type IV or type V. If A = then the model is either a Bianchi type I or a 
Bianchi type II model. 

There is one constraint equation that must also be satisfied: 



G(X) = EiV - 



A 2 - AT? + = 0, 



(2.16) 



Therefore the state space is six-dimensional; the seven evolution equations (2.5)-(2.11) are subject to the constraint 
equation ( 2.16 ) . We shall refer to the seven-dimensional state space ( [2.4| ) as the extended state space. 

By definition A is non-negative, which implies from equations ( 2.16| ) and (2.13) that E and N are also non-negative. 
Thus we have 



A > 0, 



S > 0, 



N > 0. 



(2.17) 



In addition, from the physical constraint O > together with equation (2.14), we find that the state space is compact. 
Indeed, we have that 



< \ Z 2 + ,Y,,A 2 ,A,N 



< 1. 



(2.18) 



Since both A and N are bounded, we have from equation (2.13) that N + is bounded. In equation ( pTq ) we take the 
"positive square root". In principle, there exists negative and positive values for $, but from the definition ( [2.3] ) a 
negative $ implies a negative H and hence H < for all time; i.e., the models are contracting. Since the system is 
invariant under $ — > — $, without loss of generality we shall only consider $ > 0. 



A. Invariant Sets 



There are a number of important invariant sets. Recall that the state space is constrained by equation (2.16) to be 
a six-dimensional surface in the seven-dimensional extended space. Taking the constraint equation ( 2.1 6| ) into account 



we calculate the dimension of each invariant set. These invariant sets can be classified into various classes according to 
Bianchi type and/or according to their matter content. Some invariant sets (notably the Bianchi invariant sets) have 
lower-dimensional invariant subsets. Equilibrium points and orbits occuring in each Bianchi invariant set correspond 
to cosmological models of that Bianchi type. The notation used here has been adapted from p5fl . Various lower- 
dimensional invariant sets can be constructed by taking the intersection of any Bianchi invariant set with the various 
Matter invariant sets. For example, B(I) n Ai is a 3-dimensional invariant set describing Bianchi type I models with 
a massless scalar field. 



Bianchi Type 


Notation 


Dimension 


Restrictions 


Bianchi I 


BOO 


4 


A = A = N+ = 




S(I) 


2 


A = E+ = E = A = iV + =0 


Bianchi II 


B ± (H) 


5 


,4 = 0, N+ > or N+ < 




S ± (II) 


4 


1 = 0, S = 3E+ 2 , A = E+A>+ 


Bianchi IV 


B ± (IV) 


6 


1 = 0, A>0, N+ > or N+ < 


Bianchi V 


B(V) 


4 


1 = 0, A > 0, £+ = A = 7V+ = 




S(V) 


3 


1 = 0, A>0, E + = E = A = N+ = 


Bianchi VIh 


B(VI h ) 


6 


l<0, A > 




S(VI h ) 


4 


I < 0, A > 0, 3E^ + ZE = 0, iV+ = A = 




S ± (III) 


5 


l = -l, A>0, 3T, 2 + — E = 0, A = E+JV+ 


Bianchi Vllh 


B ± (VII h ) 


6 


I > 0, A>0, N+ > or N+ < 




S ± (VII h ) 


3 


l>0, A>0, E+ = E = A = 0, Nl=lA>0 



TABLE I. Bianchi Invariant Sets. We note that B(I) and B ± (II) are class A Bianchi invariant sets which occur in the 
closure of the appropriate higher-dimensional Bianchi type B invariant set (see Fig. 1). In addition, if I is non-negative, 
N+ > and N+ < define disjoint invariant sets (indicated by a superscript ± in the table). Due to the discrete symmetry 
A — > —A, 7V+ — > — N+, these pairs of invariant sets are equivalent. 
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Matter Content Notation Dimension Restrictions 



Scalar Field 


5 


5 


Q 


= 




$ 


7^0 


Massless Scalar Field 


M 


4 


n 


= 




$ 


= 


Vacuum 


V 


3 


n 


= 


* = 0, 


$ 


= 


Perfect Fluid + Scalar Field 


TS 


6 


Q 






$ 




Perfect Fluid + Massless Scalar Field 


TM 


5 


Q 






$ 


= 


Perfect Fluid 


T 


4 


n 




* = 0, 


$ 


= 



TABLE II. Matter Invariant Sets. 



An analysis of the dynamics in the invariant sets V and T has been presented by Wainwright and Hewitt |p6| . 
Equilibrium points and orbits in the invariant set M. correspond to models with a massless scalar field; i.e., scalar 
field models with zero potential. These models are equivalent to models with a stiff perfect fluid (i.e. 7 = 2) equation 
of state; see [ p6[ . Equilibrium points and orbits in the invariant set TM. can be interpreted as representing a two- 
perfect-fluid model with 72 = 2 p7f . A partial analysis of the isotropic equilibrium points in the invariant set iS was 
completed by van den Hoogen et al. pq ]. We note that the so-called scaling solutions [^l],^9],[lo| are in the invariant 
set TS. 

The isotropic and spatially homogeneous models are found in the invariant sets S^iVIIh) U S(I) if I 7^ 0, and 
S(V) U S(I) if I = 0. In particular the zero curvature isotropic models are found in the two dimensional set S(I), 
while the negative curvature models are found in the three-dimensional sets S^(VIIh) or S(V) depending upon the 
value of I. See van den Hoogen et al. for a comprehensive analysis of the isotropic scaling models p3| . 

We note that in the invariant set B(I) there exists the invariant set S + + "f 2 < 1, A = A = N + = <t> = 0, which 
may be directly integrated to yield 

-1 

, C — constant, (2.19) 

where r is the time parameter. This solution asymptotes into the past towards the paraboloid K, (section III.B.l), and 
asymptotes to the future towards the point P(7). This solution belongs to the matter invariant set J-AA, asymptoting 
into the past towards the set M . 



1 + Ce 
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B. Monotone Functions 



The existence of strictly monotone functions, W(X.) : K" — > M, on any invariant set, S, proves the non-existence of 
periodic or recurrent orbits in S and can be used to provide information about the global behaviour of the dynamical 
system in S. (See Theorem 4.12 in [^5| for details.) 



Function: Wi(X) 


Derivative: W[ (X) 


Region of Monotonicity 


Wi = {l + Z+f-A 


Wi = -2 (2 - q) Wi + 3 (1 + E+) (2$ 2 + (2 - 7) ti) 


Monotonically approaches zero 
in the invariant set M U V. 


W 2 = 


Wi = -W 2 (2 - 3 7 ) - i (e 2 + + E) 


Monotonically decreasing to zero 

in the set (TS U TM U T)\S{I) 
when < 7 < 2/3 


Q 




Wi = -2 (2 - q) W A - 4(AiV + + E+A) 


Monotonically decreasing to zero 

in the invariant sets B(I)\S(I) and B{V)\S(V). 


A 2 


^ = 3m(, + 2 E ^- A ) 


Monotonically approaches zero 

in the invariant set S ± (I I I)\(S U J-S) , 

when 7 > 2/3. 



TABLE III. Functions, their derivatives and the sets in which they are monotonic. 
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Hewitt and Wainwright found a number of monotone functions in the invariant sets of dimension less than four in 
the perfect fluid case (i.e., in lower-dimensional subsets of the perfect fluid invariant set) and these are summarized 
in an Appendix in Hewitt and Wainwright p5 26 1. However, they were not able to find a monotonic function in the 
full perfect fluid invariant set for 2/3 < 7 < 2. 



C. The Constraint Surface 



The constraint equation G(X) = and the Implicit Function Theorem can generally be used to eliminate one of 
the variables at any point in the extended state-space provided the constraint equation is not singular there, i.e., 
grad(G(X)) 7^ 0. The constraint surface is singular for all points in the invariant sets S(I), B(V) and S(VIIh) and 
therefore cannot be used to eliminate one of the variables (and hence reduce the dimension of the dynamical system 
to six). 

Therefore, we cannot determine the local stability of equilibrium points in the sets S(I), B(V) or S{VIIh) within 
the six-dimensional state-space, and hence we are required to determine the local stability of these equilibrium points 
in the extended space, due to the singular nature of the constraint surface. This leads to further complications because 
of the limited use of the Stable Manifold Theorem. If these equilibrium points are stable in the extended state space, 
then they are stable in the six-dimensional constrained surface. However, if these equilibrium points are saddles in 
the extended state-space, then one cannot easily determine the dimension of the stable manifold within the constraint 
surface. 



III. CLASSIFICATION OF THE EQUILIBRIUM POINTS 



Let us a naly se the evolution eq uatio ns for the matter variables, namely equations (2.10) and ( [2.11 ) and the auxiliary 
equation ( 2.11 ). From equation ( 2.15 ) we find that at the equilibrium points either 



(A) n = 0, 



(3.1) 



-7 - 1. 

2 ' 



(3.2) 



In the scalar field case (^4) there is no perfect fluid present. This is the scalar field invariant set S. The equilibrium 
points and their stability will be studied in subsection III. A. These models include the massless scalar field case in 
which <f> = (V = 0), but not the vacuum case $ = \P = which will be dealt with as a subcase of the perfect fluid 
case (see below). The equilibrium p oin ts of case (^4) include the isotropic Bianchi VII^, models studied in [ 2g|. 

If, on the other hand, equation (|3.2| ) is satisfied, assuming that 7 < 2 so that q ^ 2, from equations ( 2.1C ) and 
( |2.11 ) we have that 



(Bl) * = 0,$ = 



(3.3) 



or 



(B2) * 



-V3 7 
V2k ' 



2 _ 37(2-7) 
2fc 2 ' 



(3.4) 



In case (Bl), in which both equations (3.1) and (3.2) are valid, there is no scalar field present. The perfect fluid 
subcase, which was studied by Hewitt and Wainwright [^6|, will be dealt with in subsection III.B. Note that from 
equation ( [2.11 ) $ = is an invariant set, d enot ed A4. 

The final case (B2), in which equation (3J2) is valid and neither the scalar field nor the perfect fluid is absent, 
corresponds to the scaling solutions when 7 > 0. If we define 



- 1 12 



(3.5) 



then from equation fl3.4|) we find that 



70 



2^ 2 

\[/2 1 (£2 



= 7, 



(3.6) 
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so that the scalar field "inherits" the equation of state of the fluid. It can be shown that there are exactly three 
equilibrium points corresponding to scaling solutions; the flat isotropic scaling solution described in |21||, and whose 
stability was discussed within Bianchi type VII^ models in S), and two anisotropic scaling solutions |24|. This will 
be further discussed in subsection III.C. 

Hereafter, we shall assume that < 7 < 2. The value 7 = corresponds to a cosmological constant and the model 
can be analyzed as a scalar field model with the potential V = Vo + Ke k ^ fl3l| . The value 7 = 2, corresponding to the 
stiff fluid case, is a bifurcation value and will not be considered further. 



A. Scalar Field Case 



There are seven equilibrium points and one equilibrium set in the scalar field invariant set S in which Q = 0. The 
first three equilibrium points were given in JIt]] (wherein matter terms were not included): they represent isotropic 
models (S + = E = N = A = 0): 



1) P 5 (J): E+ = E = A = i = iV + = 0,* = -fc/V6, $ = y/l - fc 2 /6 

This equilibrium point, for which q = — 1 + fc 2 /2 and which exists only for fc 2 < 6, is in the Bianchi I invariant set 
B(I). This point represents a flat FRW power-law inflationary model |q,|ri|. The corresponding eigenvalues in the 
extended state space are (throughout this paper, we shall not explicitly display the corresponding eigenvectors): 

-\{& -k 2 ), -i(6-fc 2 ), -(6-fc 2 ), -(4-fc 2 ), -(2-fc 2 ), -^(2-fc 2 ), fc 2 -3 7 . (3.7) 



2) P±(VII h ): E + = E = A = 0,i = ^1 , N + = ± vEJE 2 ! j9 = _ V| ; $ = _2_ 

These two equilibrium points (the indices "±" correspond to the ± values for N + ), which occur in the Bianchi VII/, 
invariant set S(VIIh) (since A > 0, then k 2 > 2 and therefore I > 0), have q = 0. These equilibrium points represent 
an open FRW model p8[. The corresponding eigenvalues in the extended state space are: 



2-3 7 , -1 ±^ y / k 2 - 8/3, -2±^Jk 2 - 4(fc 2 -2)l±\] [k 2 - 4(fc 2 - 2)lf + 161 (k 2 - 2) 2 + k 4 



(3.8) 



2a) P S (V): E+ = E = A = 0,i = , N + = 0, * = $ = a 



This case corresponds to points 2) for I = and belongs to the set S(V). The corresponding eigenvalues in the 
extended state space are: 

V3i 



3 7 , -1±-V- \/k 2 -8/3, -2, -2, 0, -4. (3.9) 
k 



3) P± (II): E+ = , E = 3E 2 ,, A = E + iV + , A = 0, iV + = ±3^^^, y = _J^L, $ = 6 « 

These two equilibrium points, for which g = 8(fc 2 — 2)/(fc 2 + 16) > 0, exist only for 2 < k 2 < 8. These two points 
represent Bianchi type II models analogous to those found in |2q|. The corresponding eigenvalues are: 



10 fc 2 -2 R k 2 -8 c fc 2 -8 (fc 2 -8)± V(13fc 2 -32)(fc 2 -8) fc 2 

fc 2 + 16' V + 16' V + 16' fc 2 + 16 ' 7+ V + 16' { ' 



4) P s (V7 fc ): E + = ^^,E = -3E»/i,A = 0,1= 9 ( fc2 ~y ~ 2 > ,iV + = 0, * = = W^-sOd-P , 

where n = fc 2 (Z — 3) + 41. Since E > 0, we have that I < and hence this equilibrium point occurs in the Bianchi 
Vlh invariant sets. The deceleration parameter is given by q = 2l(k 2 — 2)/[k 2 (l — 3) + 41] > 0, where k 2 > 2, and this 
point corresponds to a Collins Bianchi type VI^ solution fl32j. The corresponding eigenvalues are: 

c fc 2 ~ 2* _ . fc 2 (l-0 

DTTTTl ^ 777 > — 37 — 6- 



[fc 2 (Z - 3) + 4Z] ' ' [fc 2 (Z-3)+4Z]' 

(fc 2 - 21) ± y/(fc 2 - 2Q 2 + 81(1 - l){k 2 - 2) (fc 2 - 21) ± V(fc 2 - 2Q[(fc 2 - 2Q - 4(1 - Q(fc 2 -~2)] 

[fc 2 (Z-3) + 4Z] ' [fc 2 (Z-3)+4Z] ' [ ' 
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Let us next consider the Massless Scalar Field Invariant Set M. : there is one equilibrium set which generalizes the 
work in p(| to include scalar fields: 

5) /Cm: £ + £ + + 1' 2 = l,A = i = iV + = $ = 0,*^0 

This paraboloid, for which q = 2, generalizes the parabola K in ||(| defined by £ + £ + = 1 to include a massless 
scalar field, and represents Jacobs' Bianchi type I non- vacuum solutions p2]. However, the eigenvalues are considerably 
different from those found in [p6| , and so we list them all here (the variables which define the subspaces in which the 
corresponding eigendirections reside are included below in curly braces): 



2[(1 + E+)±V3E], 0, 0, 3(2-7), 4(1 + 53+), ^(V6 + k*). 

{A,N+} {£+,£} {£+,£,*} {£+,£,*} {£+, £, A, {£,$} 1 ' ' 



B. Perfect Fluid Case, * = $ = 

As mentioned earlier, the perfect fluid invariant set T in which if? = $ = was studied by Hewitt and Wainwright 
p6f ; hence this subsection generalizes their results by including a scalar field with an exponential potential. We shall 
use their notation to label the equilibrium points/sets. There are five such invariant points/sets. In all of these cases 
the extra two eigenvalues associated with if? and <& are (respectively) 

-|(2- 7 )<0, ^ 7 >0. (3.13) 

1) P(I): £+=£ = A = i = AT+ = # = $ = 

This equilibrium point, for which ft — 1, is a saddle for 2/3 < 7 < 2 in T (2(| (and is a sink for < 7 < 2/3), 
which corresponds to a flat FRW model. 



2) P±(II): £+ = -i(3 7 — 2), £ = 3£ + , A = E+JV+, A = 0, N+ = ±§ ^(37 - 2)(2 - 7), * = $ = 
This equilibrium point, for which f2 = ^(6 — 7), is a saddle in the perfect fluid invariant set ptj ]. 

3) P(VIfc): E+ = -i(3 7 — 2), E = -3E 2 + A A = 0, i = -^(37 - 2)(2 - 7), iV+ = * = $ = 

Since E > and A > 0, this equilibrium point occurs in the Bianchi VI/j invariant set and corresponds to the 
Collins solution jj2), where fi = |(2 - 7) + ^(37 - 2) (and therefore 2/3 < 7 < 2(-l - l)/(3 - J) and so 2 < -1). 
In p(J] this was a sink in J 7 , but is a saddle in the extended state space due to the fact that the two new eigenvalues 
have values of different sign. 

There are also two equilibrium sets, which generalize the work in |p6[ to include scalar fields: 

4) Cf: £ = -E+(l + E+), A = 0,1 = (1 + E+) 2 , 7V+ = ±VQ_ + + £+) - 3E+], * = $ = 
For this set O = 0. The local sinks in this set occur when [E6[ 

(a) Z < (Bianchi type VL h ) for -±(3 7 - 2) < £+ < Z/(3 -7) and Z > -(3 7 - 2)/(2 - 7) < 0, 

(b) Z = (Bianchi type IV) for -i(3 7 - 2) < E+ < 0, 

(c) Z > (Bianchi type VU h ) for -\(&y - 2) < E+ < 0. 
The additional two eigenvalues for the full system are: 

1-2E+, -2(1 + E+). (3.14) 

Finally, let us consider the Massless Scalar Field Invariant Set T M: 

5) K: E + E^_ = l,A = l = iV+ = $ = * = 

This parabola, for which q = 2, is the special case of ICm for which ^ = and corresponds to the parabola K. in 
p6| |. However, the eigenvalues are considerably different from those found in |2(| and so we list them all here (the 
variables define the subspaces in which the corresponding eigendirections reside are included below in curly braces): 



2[(1 + E+)±V3E], 0, 0, 3(2-7), 4(1 + E+), 3. 

{A,A+} {E + ,E} m {£+,£} {£+,£, i} {$} 1 ' > 
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We include here the equilibrium points/sets and corresponding eigenvalues as listed in |2q ]. 



Eqm. point /set 






Eigenvalues 






Comment 


P(I) 


-1(2-7) 


-3(2 


- 7) (37 - 


4) 


i(3 7 -2) i(3 7 - 




2) 




P±(77) 


1(37-2) 


-1(2- 


7) -1(2- 


-7)Wl-^ 


-7) | 


Constraint eqn. used to eliminate E 


P(V/k) t 


-1(2- 


7)(1± VI - r 2 ) - 


|(2- 7 )(i± -g 2 : 






Constraint eqn. used to eliminate E 


K. 


2(1- 


-s+) 


2(2-7) 2 


H 


-E+i^l-E^) 


1-D invariant set 


V 





2 


1 + E+ ± \/3E 


2(1 + E+) 






2-D invariant set, 7 = 2 


Ct 


- 


-4E+- 


(37 " 2) - 


- 2 


(l + E+)±2iiV + ] 






Constraint eqn. used to eliminate E 








-2 4 - 2 


3 






/ > 0, non-hyperbolic, 7 = 2/3 



TABLE IV. Equilibrium sets found by Hewitt and Wainwright [26], and the corresponding eigenvalues in the extended space. 
In the table r 2 = 2(3 7 - 2)(1 - l c /i), q 2 = 2r 2 /(2- 7) and l c = -(37 - 2)/(2 -7). 
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C. Scaling Solutions 



Defining 



2k 1 



<3> 2 



37(2 - 7) 
2fc 2 ' 



(3.16) 



and recalling that < 7 < 2, there are three equilibrium points corresponding to scaling solutions. Because the scalar 
field mimics the perfect fluid with the exact same equation of state (7^ = 7) at these equilibrium points, one may 
combine these two "fluids", via p tot = p^ + p, [i tot = \i§ + fx, p tot = (7 — l)^toti therefore, all of these equilibrium 
points will correspond to exact perfect fluid models analogous to the equilibrium points found in pq| . 

The flat isotropic FRW scaling solution ||© 
1) F S (I): E+ = E = A = A = N+ = 0, * = $ = $ s 

The eigenvalues for these points in the extended space, for which £1 = 1 — 37/fc 2 (and therefore k 2 > 37) are: 



1 



--(2-7), -3(2-7), 3 7 -4, 3 7 -2, -(3 7 - 2), --(2 - 7) ± 4^(2 - 7)(2 - 9 7 + 24 7 //c 2 ) 
There are two anisotropic scaling solutions: 



(3.17) 



2) As(II): S+ = -i(3 7 — 2),E = 3E 2 ,, A = E+JV+, A = 0, N + = ±§^(37 - 2)(2 - 7), * = * s , $ = <& s 
The eigenvalues for these points, for which = ^(6 — 7) — 3j/k 2 (and therefore fc 2 > 167/(6 — 7]), are: 



4(37-2) 



■2(2-7), 



(2-7)±W(2-7) 2 -^(2-7)^Oi: 



7(6 - 7) 
k 2 



I67 

6 — 7 



(3.18) 



where E\ 



2(37-2)- ^pl (k 2 -^L)]\ |(37-2)2^ (** " ^) ■ 



3)^ S (F4): E+ = -I (3 7 - 2) , E = -3E 2 ,//, A = 0, i = -^(2 - 7) (3 7 - 2) , iV + = 0, * = * s , $ = $ s 

These points occur in the Bianchi VI?,, invariant set (I < since E > 0) for which fl = |(2 — 7) + |j(37 — 2) — 3j/k 2 
(and therefore — l^ 1 < (2 — 7)/(37 — 2) and k 2 > 47/[(2 — 7) + (37 — 2)/l] ) and correspond to the Collins Bianchi 
Vlh perfect fluid solutions p2fl. The eigenvalues for these equilibrium points are: 



(2- 7 )± J(2- 7 ) 2 -4(3 7 -2)- 



1 



( 2-7 
\37-2 / 



where E2 = 



(2 - 7) ± y (2 - 7) 2 - * ^ - 'J 
47 (1 - |?) - (37 - 2 



2 -7 
37-2 



4:(l"|f)+(*>,-2, 



128^. 



/ 2-7 
V37-2 



(3.19) 



IV. STABILITY OF THE EQUILIBRIUM POINTS AND SOME GLOBAL RESULTS 

The stability of the equilibrium points listed in the previous section can be easily determined from the eigenvalues 
displayed. Often the stability can be determined by the eigenvalues in the extended state space, otherwise the 
constraint must be utilized to determine the stability in the six-dimensional state space (i.e., within the constraint 
surface) . In the cases in which this is not possible, we must analyse the eigenvalues in the extended seven-dimensional 
state-space, and the conclusions that can be drawn are consequently limited. Employing local stability results and 
utilizing the monotone functions found in Table Ell w e are able to prove some global results. In the absence of 
monotone functions, and in the same spirit as Refs. p5| and p6[ , we conjecture plausible results which are consistent 
with the local results and the dynamical behaviour on the boundaries and which are substantiated by numerical 
experiments. 
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A. The Case Q = 



If fl = and $ = 0, then the function W\ in Table [II monotonically approaches zero. The existence of the monotone 
function W\ implies that the global behaviour of models in the set M U V can be determined by the local behaviour 
of the equilibrium points in M. U V. Consequently, a portion of the equilibrium sets K, and ICm (corresponding to 
local sources) represent the past asymptotic states while the future asymptotic state is represented by or in the 
case of Bianchi types I and II, by a point on /C. 

Therefore, all vacuum models and all massless scalar field models are asymptotic to the past to a Kasner state and 
are asymptotic to the future either to a plane wave solution (Bianchi types IV, VI/j and VII/j), or to a Kasner state 
(Bianchi types I and II), or to a Milne state (Bianchi type V). 

If ft = and $ ^ 0, then the models only contain a scalar field. It was proven in (ll| that all Bianchi models 
evolve to a power-law inflationary state (represented by Ps(J)) when k 2 < 2. If k 2 > 2, then it was shown in pSy 
that a subset of Bianchi models of types V and VII^ evolve towards negatively curved isotropic models represented 
by points P S (V) and P±{VII h ). In || it was shown that when k 2 > 2 the future state of a subset of Bianchi type 
Vlh solutions is represented by the point Ps(VIh)- It can be seen here that the future state of a subset of Bianchi 
type II models is represented by the point Pg(II). 

Therefore, all scalar field models with f2 = evolve to a power-law inflationary state if k 2 < 2. If k 2 > 2, then the 
future asymptotic state for all Bianchi type IV, V and VII/j models is conjectured to be a negatively-curved, isotropic 
model and the future asymptotic state for all Bianchi type Vlft, models is conjectured to be the Feinstein-Ibanez 
anisotropic scalar field model [fL4[| . If 2 < k 2 < 8, then the future asymptotic state for all Bianchi type II models is 
the anisotropic Bianchi type II scalar field model, and if k 2 > 8 then the future asymptotic state is that of a Kasner 
model. If 2 < k 2 < 6, then the Bianchi type I models approach a non-inflationary, isotropic (i.e., the point Ps(/)); if 
k 2 > 6, then they evolve to a Kasner state in the future. 

B. The Case Q / 0, < 7 < 2/3 



If O 5^ and < 7 < 2/3 then the function W2 in Table III is monotonically decreasing to zero. Therefore, 



we conclude that the omega-limit set of all non- exce p tiona l orbits (i.e., those orbits excluding equilibrium points, 



heteroclinic orbits, etc.) of the dynamical system (2.5)-(2.11) is a subset of S(I). This implies that all non-exceptional 



models with Q =/= evolve towards the zero-curvature spatially-homogeneous and isotropic models in S(I) and hence 
isotropize to the future. In ]23( ], it was shown that the zero-curvature spatially-homogeneous and isotropic models 
evolve towards the power-law inflationary model, represented by the point Ps(-0 when k 2 < 37 or towards the isotropic 
scaling solution, represented by the point Fs{I), when k 2 > 37. Using Wi, we also conclude that the past asymptotic 
state(s) of all non-exceptional models (including models in S(I)) is characterized by f2 = 0. In other words, matter is 
dynamically unimportant as these models evolve to the past. It was shown in p3] that all models evolve in the past 
to some portion of K, or K-m (the Kasner models) which are local sources. 

C. The case fi / 0, |<7<2 
The following table lists the local sinks for | < 7 < 2. 
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Sink 


Bianchi type 


k 


Other constraints 


p s (vn h y 


I 
I 


k 2 < 2 
k 2 =2 




Ps(VI-i) 


III 
III 


k 2 >2 
all 


7 > fcV^ + 1), Z = -l 
7 > 1, S+ = -1/4 


?s(V) 


V 


k 2 > 2 




Ps(VI h ) 
As{VI h ) 


VI h 
Vlh 
VI h 


k 2 > 2 
all 


j>2k A (l-l)/\k 2 (l-3)+4l} 

1 \ / / L \ / 1 J 

7>4/3, E+<-l/2 
7 ^ -(37-2) 


- [(2- 7 )+(3 7 -2)/il 






Vllh 


k 2 > 2 


' > 4(^-2^4-^) for 2 < fc 2 < 4 

' < 4(^-2^-4) fOT fc2 > 4 



TABLE V. This table lists all of the sinks in the various Bianchi invariant sets for 2/3 < 7 < 2. A subset of JCm acts as a 
source for all Bianchi class B models. * Note: in this case N + = (i.e., Ps=P$ ) and in fact corresponds to a Bianchi I model. 
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The function W3 is monotonically decreasing to zero in B(I)\S(I) and B(V)\S(V). This implies that there do not 
exist any periodic or recurrent orbits in these sets and, furthermore, the global behaviour of the Bianchi type I and V 
models can be determined from the local behaviour of the equilibrium points in these sets. We conjecture that there 
do not exist any periodic or recurrent orbits in the entire phase space for 7 > 2/3, whence it follows that all global 
behaviour can be determined from Table [y| 

We note that a subset of Km with (1 + S+) 2 > 3S, ^ > —y/6/k acts as a source for all Bianchi class B models. For 
k 2 < 2, Ps(-T) is the global attractor (sink). From Table |y| we see that there are unique global attractors (both past 
and future) in all invariant sets and hence the asymptotic properties are simple to determine. The sinks and sources 
for a particular Bianchi invariant set, which may appear in that invariant set or on the boundary corresponding to a 
(lower-dimensional) specialization of that Bianchi type, can be easily determined from Table M and Figure 1 which 
lists the specializations of the Bianchi class B models |33| . 



(A) 


B 


V \ / ll0> 
II ^ ' 

I 


„VL „-VII 

^^^^^/^^ 
- ' ' IV 

, V 



Figure 1: Specialization diagram for Bianchi class B models obtained by letting a non-zero parameter go to zero. A 

broken arrow indicates the group class changes (from B to A). 



The most general models are those of Bianchi types VI/j and VII?,,. The Bianchi type VII/,, models are of particular 
physical interest since they contain open FRW models as special cases. From Table |y| and Figure 1 we argue that 
generically these models (with a scalar field) isotropize to the future, a result which is of great significance. The 
Bianchi type VI/j models are also of interest since they contain a class of anisotropic scaling solutions that act as 
attractors for an open set of Bianchi type B models. We note that generically Bianchi type VI^ models do not 
isotropize for k 2 > 2. 



V. INTERMEDIATE BEHAVIOUR AND THE INVARIANT SET S(VI H ) 

It is also of interest to determine the intermediate behaviour of the models. In order to do this, we need to investigate 
the saddles, determine the dimension of their stable submanifolds, and construct possible heteroclinic sequences. This 
could then be used, in conjunction with numerical work, to establish the physical properties of the models. For 
example, we could investigate whether intermediate isotropization can occur in Bianchi type VII^ models |}4| . There 
are many different cases to consider depending upon the various bifurcation values and the particular Bianchi invariant 
set under investigation. As an example, we shall study the heteroclinic sequences in the four-dimensional invariant set 
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S(VIh), because it illustrates the method and because such a study emphasizes the importance of anisotropic scaling 
solutions. 

The subspace S(VIh), which arises from the restrictions N + = A = and 3S+ + IT, = |2j| is, in fact, the class of 
diagonal Bianchi VI^ models and is four-dimensional (and was shown in J2~if| to illustrate th e exi st ence and importance 
of the anisotropic scaling solutions). From the above restrictions, the system of equations (2.5)-( 2.11 ) now reduce to: 



§L4 



£'+ = (q - 2)E H 
A' = 2(g + 2S + )A 
= (g-2)#- \VQk<$> 2 
= (q + 1 + iy^A;*)* 



(5.1) 
(5.2) 
(5.3) 
(5.4) 



We note that the constraint equation ( 2.16| ) is automatically satisfied. 

The equilibrium points of the system are a subset of the ones presented in section III and we present those that 
belong to this subspace and their corresponding eigenvalues in Table VI: 



Eqm point/set 


Eigenvalues 


Stability 


Conditions 




6 6 oi-1 3 (3 7 -2)+/(2- 7 ) 
1-3 1-3 1-3 2 1-3 


s 2 < 37 - 2 : W s 
s 2 > 3 7 - 2 : W s 


{$ = 0} 
{$ = Q = 0} 






3(2- 7 ) 4(1 + E+) 


k 2 < 6 W s 
fc 2 > 6 and S§ > ^ IF S 
fc 2 > 6 and < £ VP 


{0} 
{0} 

{fi = A' = 0} 




P(I) 


-f(2- 7 ) -1(2-7) 3 7 -2 | 7 


IF S 


{$ = A = 0} 




P(VI h ) 


| 7 -|(2-7) -!(2- 7 )(l±Vl-r 2 ) 




{$ = 0} 


s 2 > 3 7 - 2 


Ps(I) 


fc V 6 fe ^ 6 k 2 2 k 2 3 7 


fc 2 < 2 iy s 

2 < fc 2 < 3 7 F^ s 
fc 2 > 3 7 W s 


{S(V/ h )} 
{K = 0} 

{k = n = 0} 


fc 2 < 6 


Ps(VI h ) 


[**(J-3)+4Z] ^ D [/c 2 (i-3)+4i] 
g (fc 2 -2i)±V 'O 2 -2i)[(fc 2 -2i)-4(l-()(fc 2 -2)] 

** [fe 2 (;-3)+4;l 


k 2 < 37 or s 2 < Sq W s 
s 2 > si W s 


{S(VA)} 
{fi = 0} 


fc 2 > 2 




-|(2- 7 ) 3 7 -2 
-|(2 - 7) ± V(2 - 7)(2 - 9 7 + 24 7 2 /fc 2 ) 


w s 


{A = 0} 


fc 2 > 3 7 


A s {VI h ) 


-|(2 - 7) ± Va± Vb 


W s : {S(VI h )} 


s 2 > s 2 



Table VI: This table lists all equilibrium points of the system ( 5.1)-&5 
(sub) manifold, W 



and gives the corresponding stable 



We also give the conditions on the parameters for their existence. We note that s 2 = 
S 2 = 1 - t - E 2 . = 1 - ME* r 2 = 2(3 7 - 2)(1 - l c /l), q 2 = 2r 2 /(2 - 7) and l c = -(3 7 - 2)/(2 - 7). 



2 = fc 2 (3 7 -2) 
*0 — (fc 2 -3 7 ) ' 

Note that the four eigenvalues for As(VIh) are the same as the last four eigenvalues in Eq. (3.H ), but are listed here 
in "short form" for compact notation. 



^From the eigenvalues in Table V we can study the stability of the equilibrium points and the qualitative behaviour 
in this four-dimensional subspace. The bifurcation values leading to a change in the asymptotic behaviour are: 
k 2 = 2, fc 2 = 6, /c 2 = 3 7 , s 2 = 3 7 -2, s 2 = fc 2 (37-2)/(fc 2 -3 7 ). 

We can see that the possible future asymptotic behaviour is given by points containing a scalar field, either alone 
or together with a perfect fluid component in the case of the anisotropic scaling solution. These possible future 
asymptotics are: Ps(I), Ps(VIh) and AsiVIh). 



A. Heteroclinic Sequences 



The stable and unstable manifolds of the saddle equilibrium points provide a skeleton of special orbits that play 
a significant role in determining the dynamics of the models (and in particular, their intermediate states) ]34| ]. A 
(finite) heteroclinic sequence is a set of equilibrium points E Ql E\, . . . ,E n , where E is a local source, E n is a local 
sink and the rest are saddles, such that there is a heteroclinic orbit which joins E^\ to Ei for each i = 1, ... ,71 p5[ . 
Below are the heteroclinic orbits in the class of S(VIh) models for different parameter ranges. Equilibrium points in 
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parentheses stand for optional intermediate points. We recall that s 2 = 
the Kasner-like ring splits into Km [A] (sources) and Km [B] (saddles.) 



k 2 < 2 < 3 7 



i) s 2 < 37 - 2 : 

K M -> Cf -> P 5 (J) 

/C M ^ P(7) ^ £± ^ P 5 (7) 



jIL aT1 d s 2 = fc (37-2) 
Z-3 anQ S - (fc2-3 7 ) • 



Finally, when fc 2 > 6 



ii) 37 - 2 < s 2 

K M -> £± -> P(F4) -> P s (/) 
/C^ -> -> P(^4) -> Ps(-0 



2 < fc 2 < 37 



i) 5 2 < 3 7 - 2 : 

/Cm - P(/) - Ps(7) -> P|(^4) 

/Cm-p(/)-^-pJ=(^4) 



ii) 37 - 2 < a 2 : 

(P(VI h )) -> P|(^/4) 
/Cm - P(^4) - P|(F/4) 

/Cm - P(/) - PsW - P|(^4) 



i) a 2 < 37 - 2 < s 2 , : 

K M - PsU) - (-F 5 (/)) - P|(^/4) 
/Cm-(P(/))-^-P 5 ± (^/4) 
/Cm -> P(/) - -» P|(VJ/ fc ) 

iii) 37 - 2 < a 2 < s 2 : 
/Cm -> ^fe -> P|(VJ/ fc ) -» A s (VI h ) 
K M ^£±^ P(VI h ) - 4s(W fc ) 
/Cm - Ps(/) - Pj=(VJ/ fc ) - As(^) 
/Cm - PsOO - FsV) - ./4s(W h ) 
/Cm - P(/) - P(W fc ) - ^s(VJ h ) 
/Cm - P(/) - Fs{I) -» ^5(^4) 



ii) 37 - 2 < s 2 < sg : 
K M ^Ct^ (P{VI h )) - P|(^/4) 
/Cm - Ps(/) - (Jvs(/)) - P|(VJJ fc ) 
/Cm -> P(/) -> P(W fc ) -» P|(^J/ fc ) 
/Cm - -> - Ps(VJJ fc ) 



37 < 6 < fc 2 



i) s 2 < 37 - 2 < s 



/Cm [A] 
/Cm [A] 



(/Cm[P]) 
(/Cm[P]) 



ii) 37 - 2 < s 2 < s 2 : 
K M [A]^{K M [B])- 
K M [A] - (/Cm[P]) - 
/Cm[A]^(/Cm[P])- 

iii) 37 - 2 < s 2 < s 2 

/Cm[A]^(/Cm[5])- 
/CmM^(/Cm[P])- 
/Cm[A]^(/Cm[S])- 
/Cm[A]^(/Cm[S])- 



(P(/))-£±-P|(F/4) 
P(/)^jr s (7 ) ^p±(y77 ft ) 



P(7) 



(P(^)) 
-> P(VI h ) 



p£{vii h ) 
p£{vih) 



P{I)^Ts{I)^Pf{VII h ) 



Pf{VII h )^A s {VI h ) 
P{VI h )^A s {VI h ) 
P(I)^P(VI h )^A s (VI h ) 
P{I)^ T S {I)^ A s (VI h ) 
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K M >- P S (I) P(VI h ) ► A s {VI h ) 




P{I) *■ Fs{I) 



Figure 2. Heteroclinic sequences in S(VIh) for s < 3j < k 2 < 6 and 37 — 2 < s§ < s 2 , indicating the skeleton of orbits 
defined by the stable and unstable manifolds of the saddle points. Note that the anisotropic scaling solution As{VIh) 
is a stable attractor. 

Figure 2 can be used to construct sequences of orbits joining equilibrium points starting at the past attractor K,m 
and ending at the future attractor As(VIh)- For each sequence there will be a family of orbits that shadow this 
sequence in the state space. 

VI. CONCLUSION 

In this paper we have discussed the qualitative properties of Bianchi type B cosmological models containing a 
barotropic fluid and a scalar field with an exponential potential. The most general models are those of type VI^, 
which include the anisotropic scaling solutions, and those of type VII^, which include the open FRW models. 

In cases in which we have been able to find monotone functions we have been able to prove global results. Otherwise, 
based on the local analysis of the stability of equilibrium points and the dynamics on the boundaries of the appropriate 
state space, we have presented plausible global results (this is similar to the analysis of perfect fluid models in )2if] 
and |^6| in which no monotone functions were found in the Bianchi type VI and VII invariant sets). In all cases, 
however, our results are further justified by numerical experimentation. 

Let us summarize the main results: 

• All models with k 2 < 2 asymptote toward the flat FRW power-law inflationary model |5|,[l7]], corresponding to 
the global attractor Ps(I), at late times; i.e., all such models isotropize and inflate to the future. 

• Fs(I) is a saddle and hence the flat FRW scaling solutions |^,^0| do not act as late-time attractors in general 
& 

• A subset of K-m acts as a source for all Bianchi type B models; hence all models are asymptotic in the past to 
a massless scalar field analogue of the Jacobs anisotropic Bianchi I solutions. 

• For k 2 > 2, Bianchi type VII^ models generically asymptote towards an open FRW scalar field model, represented 
by one of the local sinks Ps(V) or Pg(VIIh), and hence isotropize to the future. 

• For k 2 > 2, Bianchi type Vlh models generically asymptote towards either an anisotropic scalar field analogue of 
the Collins solution |32fl , an anisotropic vacuum solution (with no scalar field) or an anisotropic scaling solution 
p4j , corresponding to the local sinks Pg(VIh), Ck(VIh) or As{VIh), respectively (depending on the values of 
a given model's parameters - see table IV for details). These models do not generally isotropize. 

• In particular, the equilibrium point As{VIh) is a local attractor in the Bianchi VI/j invariant set and hence there 
is an open set of Bianchi type B models containing a perfect fluid and a scalar field with exponential potential 
which asymptote toward a corresponding anisotropic scaling solution at late times. 
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We should stress that our analysis and results are applicable to a variety of other cosmological models in, for 
example, scalar-tensor theories of gravity (which are formally equivalent to general relativity containing a scalar field 
with an exponential potential) M^H-B^I, theories with multiple scalar fields with exponential potentials Q and 
string theory dlH ■ 

In future work we shall study spatially homogeneous models with positive spatial curvature J4(| and Bianchi models 
of type A |fjl[ . However, our ultimate goal is to extend the techniques used in this paper and study the more interesting 
(and physically more relevant) case of spatially inhomogeneous models. 
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